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O , We present results of local spin density approximation (LSDA) pseudopotential calculations for 

the ferroelectromagnet, yttrium manganite (YMnOa). The origin of the differences between fer- 
roelectric and non-ferroelectric perovskite manganites is determined by comparing the calculated 
properties of yttrium manganite in its ferroelectric hexagonal and non-ferroelectric orthorhombic 
phases. In addition, orthorhombic YMnOa is compared with the prototypical non-ferroelectric man- 
ganite, lanthanum manganite. We show that, while the octahedral crystal field splitting of the cubic 

f~| , perovskite structure causes a centro-symmetric Jahn- Teller distortion around the Mn''"'" ion, the 
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markedly different splitting in hexagonal perovskites creates an electronic configuration consistent 
with ferroelectric distortion. We explain the nature of the distortion, and show that a local magnetic 
moment on the Mn^"^ ion is a requirement for it to occur. 
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I. INTRODUCTION 



^ ■ The observation of colossal magnetoresistance (CMR) in (La,Ca)Mn03 ]l|] has prompted a flurry of recent research 

7^ on this material and related perovskite structure manganites i^. The majority of the recent research has focused 
^ on manganites in which the large (A-site) cation is a rare earth from the left hand side of the lanthanide series. 
S^ The manganites of these large rare earth ions (lanthanum manganite through dysprosium manganite) all crystallize 

d in the cubic perovskite structure (Figure |l|), with the same low temperature orthorhombic distortion and A- type 
antiferromagnetic ordering of the Mn'^^ ions. 

During the recent studies, many rare earth manganites have been found to show strong coupling between their 
magnetic and structural, or magnetic and electric order parameters. For example, a magnetically induced structural 
phase transition has been observed in Lao.ssSro.iyMnOs P| indicating strong coupling between the local magnetic 
spin moments and the lattice structure. In Ndo.sSro.sMnOs, strong coupling between the magnetic spin moments 



and the electronic charge carriers was demonstrated when an electronic metal-insulator transition was induced by an 
external magnetic field y. Indeed the large change in conductivity with applied magnetic field, which causes colossal 
magnetoresistance is believed to originate from a similar type of phase transition B . 

The strong coupling between magnetic and electric order parameters is of particular interest in the manganites of 
yttrium and the smaller (holmium through lutetium) rare earths at the right hand side of the lanthanide series. In 
these materials, the hexagonal perovskite structure (Figure 0) is the lowest energy structure. However the metastable 
orthorhombic cubic perovskite phase can be regained by annealing at high pressure [|6|, or by appropriate choice of 
synthetic method. The phase transition between the stable hexagonal phase and the metastable orthorhombic phase 
is reconstructive and first order. 

In the hexagonal manganites, as in the cubic manganites, superexchange between adjacent Mn^+ ions causes anti- 
ferromagnetic ordering. However in addition to the magnetic ordering, the hexagonal manganites also undergo a phase 
transition to a non-centrosymmetric low temperature state which has ferroelectric ordering. Therefore the hexagonal 
yttrium and rare-earth manganites belong to the rather small group of materials known as ferroelectromagnets, in 
which magnetic and ferroelectric ordering occurs simultaneously. Although there are thirteen point groups which 
permit the occurrence of spontaneous magnetization and spontaneous polarization in the same phase pi , remarkably 
few ferroelectric materials show any kind (ferro- ferri- antiferro-, etc.) of magnetic ordering. Understanding the fun- 
damental physics that tends to discourage simultaneous ferroelectricity and magnetic ordering is an area of current 
research activity p|. And the fact that the ferroelectric hexagonal manganites can be readily transformed to their 
non-ferroelectric cubic form makes them ideal prototypes for studying this phenomenon. 

The goal of the study described in this paper is to explore theoretically the similarities and differences between 
non- ferroelectric cubic manganites and their ferroelectric hexagonal counterparts. To achieve this goal we evaluate the 
electronic and magnetic properties of some representative perovskite manganites using a plane wave pseudopotential 
(PWPP) implementation of density functional theory (DFT) within the local spin density approximation (LSDA). 
Our choice of theoretical approach is influenced by successful first-principles studies of non-magnetic perovskite 
ferroelectrics [^|0| , and by earlier work which elucidated the origins of the differences between related cubic perovskite 
manganites |l^Jl^. We calculate the differences between the electronic properties of cubic and hexagonal YMnOs to 
determine the reasons for the lack of ferroelectricity in cubic manganites, and its presence in hexagonal manganites. 
We also compare our calculated results for cubic YMnOa with those for the prototypical cubic perovskite manganite, 
LaMnOa, to establish whether the metastability of the cubic phase of YMnOs is merely the result of the small size of 
the yttrium ion, or if electronic factors also play a role. 

The remainder of this paper is organized as follows: In Section n we summarize the experimental results for 



YMnOa from the literature. In Section |II| we briefly describe our plane wave pseudopotential implementation of 

our 



density functional theory and outline the technical details of the calculations used in this work. In Section [V 



results for cubic and hexagonal YMnOa are presented, followed by a comparison of the properties of cubic YMnOa 
with those of LaMnOs. Our results are summarized in Section ^ 



II. SUMMARY OF EXPERIMENTAL RESULTS 

The literature on the hexagonal manganites is much more limited than that on the cubic perovskite manganites, 
and modern work focuses largely on yttrium manganite, YMnOs. Although yttrium is not a rare earth, it forms 
a stable trivalent cation with a similar ionic radius to those of the smaller rare earth ions. It is favored over the 
smaller rare earth ions in research studies for two reasons. First, study of magnetic ordering on the Mn'^+ ions is 
more straightforward in YMnOa than in the rare earth manganites since there is no perturbation from the rare earth 
/-electron magnetic moment. In addition, YMnOs forms both the hexagonal and orthorhombic phases on solution 
synthesis at ambient pressure |13[| making it particularly appealing for use in a comparative study. 

Early (1960s) work on YMnOs p3[ established the hexagonal phase to be ferroelectric, with the hexagonal perovskite 
structure and the PQ^cm space group. The hexagonal perovskite structure consists of ABCACB... stacking of close 
packed O layers, with the Mn ions occupying 5-fold coordinated sites, and the rare earth atoms in 7-fold coordinated 
interstices. The ideal high temperature paraelectric structure consists of two formula units per unit cell and is shown 
in Figure |[ Since this structure does not occur in the perovskite ferrites it has been suggested that the ability of 
Mn'^+ to form 5-fold trigonal bipyramids through dsp^ hybrid bonding might be a requirement p3| . In addition (not 
shown in the Figure) at low temperature the MnOs bipyramids are slightly rotated around the axis passing through 
the Mn and parallel to one of the triangular base sides. Early reports of a weak parasitic ferromagnetism iQ were 
soon shown to be the result of Mn304 impurity contaminating the powdered phase |15| , and weak ferromagnetism 
was shown to be absent in single crystal samples |Iq] . 

Modern studies of hexagonal YMnOa have established that the ferroelectricity occurs along the c axis (i.e. [0001]) 
|Qq,0, with a spontaneous P '^ 5.5 /iC/cm^ In addition a strong coupling between the ferroelectric and magnetic 
ordering has been revealed [|l8[. Although the ferroelectric Curie temperature (914K) is quite different from the Neel 
temperature (80K), anomalies in the dielectric constant and loss tangent near the Neel temperature are indicative of 
coupling between the ferroelectric and antiferromagnetic orders. The coupling has been attributed to changes in the 
phonon spectrum associated with the antiferromagnetic transition pq | . Raman and infra-red spectroscopy of the high 
temperature paraelectric and low temperature ferroelectric phases show only weak bands in the ferroelectric phase due 
to the non-centrosymmetricity. This indicates that the structural differences between the ferroelectric and paraelectric 
phases are very small. Non-linear optical spectroscopy M] shows two types of optical second harmonic spectra of 
the Mn'^^ ions in hexagonal YMnOa, one from the non-centrosymmetric ferroelectric ordering of charges, and the 
other caused by the centrosymmetric antiferromagnetic ordering of spins. Partial overlapping between the electronic 
transitions results in a nonlinear optical polarization which depends on two order parameters. A generalized Ginzburg- 
Landau formalism has been developed which shows that the second harmonic generation susceptibility is directly 
proportional to the bilinear combination of both order parameters [E0[. Finally, epitaxial thin films of hexagonal 
YMnOa have been grown on (111) MgO and (0001) ZnO:Al/(0001) |l6|, and explored for use as non-volatile memory 
devices. YMnOs presents some technological advantages over many common ferroelectric perovskites, including lower 
dielectric constant (~ 20 at room temperature) and non-volatile constituent elements. 



In the 1970s, work began on the orthorhombic phase of YMnOa. Susceptibihty measurements estabhshed that 
orthorhombic YMnOs is an antiferromagnct with a Neel temperature of 42K, although significant deviations from 
Curie- Weiss behavior were seen weh above the manganese ordering temperature M. A more detailed study pl| ] 
concluded that the Mn'^+ ordering is helical with a propagation wavevector k = [0,0.0786,0] and a helical angle from 
plane to plane of 14°. Within each plane the ordering is antiferromagnetic. 

There has also been some modern work on the orthorhombic phase of YMnOa . Attempts have been made to find a 
synthesis which will produce a larger proportion of the metastable orthorhombic phase ||22| . In addition the Raman 
spectrum of orthorhombic YMnOs has been measured and the Raman active phonon modes determined and compared 
with those of LaMnOs Q. 

III. COMPUTATIONAL TECHNIQUE 

The calculations described in this work were performed using plane wave pseudopotential (PWPP) implementations 
pjl of density functional theory pa| within the local spin density approximation (LSD A). The accuracy and efhciency 
of ab initio pseudopotential calculations (compared with all-electron calculations) is now well established for spin- 
polarized magnetic systems, and the PWPP method has been applied to a wide range of magnetic materials Pq] . 

In this work we use both the optimized pseudopotentials developed by Rappe et al. E^ and the Vanderbilt ultra-soft 
pseudopotentials [ p8[ , to reduce the energy cut-off of the plane wave expansion to around 60 and 35 Ry respectively. 
For La and Y, we constructed scalar-relativistic pseudopotentials, and for Mn and O non-relativistic pseudopotentials 
were used with the partial non-linear core correction scheme of Louie et al. p9| applied to the Mn ions. The 
pseudopotentials were tested for transferability by comparing with all-electron calculations for a range of typical 
atomic and ionic configurations. The pseudo-eigenvalues and total energies were found to be equivalent to the all- 
electron values to within a few meVs. All pseudopotentials were converted into the usual Kleinman-Bylander form 
pO{ , and the absence of ghost states was confirmed using the ghost theorem of Gonze, Kackell and Scheffler |^1|]. The 
pseudopotential construction and the total energy calculations used the Perdew-Zunger parameterization |32| of the 
Ceperley- Alder exchange correlation potential 133] with the von Barth-Hedin interpolation formula |3J]. 

The total energies and band structures were calculated using the CASTEP 2.1 P5| , p6| and PWSCF ultra-soft 
pseudopotential programs. For cubic perovskites we used a plane wave cut off of 60 Ry, which corresponds to around 
3500 plane waves in a cubic unit cell with lattice constant of around 4 A . A 6x6x6 Monkhorst-Pack |3^] grid was used 
for all calculations for cubic systems. This led to 10 k-points in the irreducible Brillouin Zone for the high symmetry 
structures, and a correspondingly higher number for the structures with lower symmetry. In calculations for the 
hexagonal structure, a 5x5x4 Monkhorst-Pack grid was used, resulting in 18 k-points in the irreducible Brillouin Zone 
for the high symmetry structures. A variable Gaussian broadening between 1 eV and 0.002 eV was applied to the 
k-point sampling to speed convergence for metallic systems. For density of states (DOS) calculations we calculated 
first-principles eigenvalues of a large (« 100 in the irreducible Brillouin zone) k-point set, then interpolated using the 
interpolation scheme of Monkhorst and Pack |3^ . We then applied the Gilat-Raubenheimer method [Bq] to integrate 



over this fine mesh. Finally, for the band structure plots, and for use in the tight-binding analysis, symmetry labels 
along high-symmetry lines were assigned using projection operators for the corresponding irreducible representations. 

IV. RESULTS 

A. Calculated Electronic Properties of cubic and hexagonal YMnOs 

In this section we compare the electronic properties of YMnOa in its metastable cubic perovskite structure with 
those of the stable hexagonal phase. We begin by calculating the electronic structure for the ideal high symmetry 
structures, without including magnetic effects, then lower the magnetic symmetry to the hypothetical ferromagnetic 
and observed antiferromagnetic phases. Finally we introduce structural distortions to simulate the effects of the onset 
of ferroelectricity. This ability to isolate structural and magnetic distortions is unique to computational studies, and 
allows identification of the essential microscopic interactions which cause the observed macroscopic behavior. 

1. High Symmetry Param,agnetic Structures 

First we present the calculated electronic structures for YMnOs in the ideal cubic and hexagonal structures, without 
spin polarization (we call this the paramagnetic (PM) phase). 

a. Cubic Figure ^ shows the calculated density of states (DOS) for cubic paramagnetic YMnOs. We chose a lattice 
constant of 3.84 A which corresponds to the average of the lattice constants in the low temperature orthorhombic unit 
cell. The plotted energy range is from -8 eV to 4 eV, and the lower lying semi-core states have been omitted for clarity. 
The Fermi level is set to zero. The broad series of bands between approximately —2 and —7 eV arises primarily from 
the oxygen 2p orbitals. Above the oxygen 2p bands, and separated from them by an energy gap, are the Mn 3d bands. 
The Mn 3d bands are divided into two sub-bands - the lower energy i2g bands, and the higher energy Cg bands - as 
a result of the crystal field splitting by the octahedral oxygen anions. The Fermi level lies near the top of the Mn 
3d t2g bands and is in a region of high density of states. The large DOS at the Fermi level confirms that the cubic 
PM structure is unstable, and that a lower energy structure could be achieved by allowing spin-polarization and/or 
structural distortion. 

b. Hexagonal The hexagonal paraelectric structure has three lattice parameters; the usual a and c parameters 
plus an internal u that gives the distance (in units of c) between oxygen and yttrium layers, and is not fixed by the 
symmetry. In ideal hexagonal structures (without tilting of the oxygen octahedra) the atomic positions in crystal 
coordinates are: Y at (0,0,0) and (0,0, i), Mn at (|,f ,3) and (|,5,|) O at (0,0,^), (0,0, |), (i,|,u),(|,i,l - u), 
(|,i,i-|-u) and (i,|,i — u), with u not fixed by the symmetry. By energy minimization we obtain a=3.518 A, c=11.29 
A, and M=0.084, in very good agreement with the experimental values a=3.539 A, c=11.3(4) A, and u=0.084. This 
suggests that the oxygen tilting which we are neglecting in these calculations indeed has little effect on the structural 
properties. Note however that in the ideal, non-magnetic, high symmetry structures, we find that the hexagonal phase 
is not energetically stabilized over the cubic phase. This indicates that the low temperature structural distortions 



and magnetic polarizations which we investigate in the next section are responsible for the stability of the hexagonal 
phase. 

In Figure we show the calculated orbital-resolved DOSs for the non-magnetic hexagonal phase at the LSDA 
minimum energy lattice parameters. For clarity only the DOS on the transverse oxygens is shown since these are 
most strongly involved in chemical bonding with the Mn ion and will be most important in our later analysis of 
ferroelectricity. As in the cubic case, the non-magnetic hexagonal phase shows a large DOS at Ep, mostly due to the 
Mn d states located in a narrow region around E^?. In sharp contrast with the cubic phase, however, the hexagonal 
crystal field splits the d states into three groups. The dxy and d^^^y^ orbitals lie on the (0001) plane and are nearly 
degenerate (only one of them is shown in the Figure). These orbitals point towards the nearby in-plane oxygens (Op), 
and so are able to participate in covalent bonding, therefore they form a rather broad band. The d^z and dyz orbitals 
are also almost degenerate. Their DOSs overlap with those of d^y and d^2_y2, but are localized in a narrower energy 
region, since they do not point directly towards oxygens and therefore their bands are not broadened by hybridization. 
Each orbital of the two doublets is roughly half occupied. Finally, dz2 (shaded gray in the Figure for clarity) is the 
highest in energy and the least occupied. It points towards the transverse oxygens (Ot) which are the closest ligands 
to the Mn ions and therefore cause the strongest crystal field destabilization. (The distance between the Mn ion and 
the transverse oxygens along the c axis is 1.875 A and that between the Mn ion and the in-plane oxygens is 2.03 A 
.) However it is important to note that a certain amount of d^'^ charge is hybridized with the Ot Pz states located ~ 
5 eV below Ep. 

In the next section we will show that these important differences in the splitting of the Mn d bands resulting from 
the different Mn coordination in the two structures, which are apparent already in the high symmetry paramagnetic 
phases, have a profound effect on their subsequent electronic, structural and magnetic properties in the realistic 
spin-polarized structures. 

2. High symmetry structural phases with magnetic polarization 

Next we calculate the properties of both cubic and hexagonal phases, still with the ideal structural high symmetry, 
but this time in the experimentally observed antiferromagnetic ground state. In both cases, allowing spin polarization 
lowers the total energy, consistent with the Stoner model, and introduces magnetic moments onto the Mn ions. 

a. Cubic In Figure ra we show the orbital resolved densities of states on single Mn and O ions, in ideal cubic 
A-type antiferromagnetic YMnOa at the same lattice constant used for the paramagnetic calculations. The A-type 
AFM phase consists of ferromagnetically aligned (100) planes of Mn ions coupled antiferromagnetically to each other. 
Due to the exact AFM symmetry, each atom has a corresponding atom in the cell with anti-aligned spin, so the up and 
down components of DOS are exchanged. The upper plot shows the up-spin and down-spin Mn 3d densities of states 
above and below the x-axis respectively, and the lower plot shows the up-spin oxygen density of states. The down-spin 
oxygen density of states is not shown since it is very similar to the up-spin. The oxygen density of states in this 
region derives almost entirely from the oxygen 2p orbitals; there is negligible contribution from the O 2s orbitals in 



the region shown. The densities of states of the dxy, dxz and dyz (not shown) orbitals are very similar. In the absence 
of the antiferromagnetic ordering these would be identical by symmetry and would form the t2g manifold. Similarly 
the dz2 and dx2_y2 bands, from the Cg manifold, are alike. The t2g band is rather narrow for both spin directions, and 
clearly shows an exchange splitting of around 2 eV. This splitting causes the majority i2g band to be fully occupied 
while the minority t2g band contains only a small amount of electron density. The Cg band is much broader since 
the 6^22 and (i^2_j,2 orbitals are more able to hybridize with their oxygen neighbors. It is roughly half-occupied in the 
majority case, but completely unoccupied in the minority case. The calculated Mn magnetic moment of around 3 
/is, significantly reduced from the Hund's rule value of 4 ^b expected for a d'^ ion, reflects the fact that the exchange 
splitting is not large enough to produce full spin polarization. It is interesting to note that the region of overlapping 
energy between the Mn 3d and oxygen 2p states is in fact rather small, with only a small amount of Mn density 
occurring in the oxygen rich region between around -3 and -9 eV, and vice versa. The similarity between the up- and 
down- spin oxygen densities of states is a a consequence of the small Mn 3d - O 2p overlap. 

For comparison, the orbital resolved densities of states in ideal cubic ferromagnetic YMn03 are shown in Figure ra. 
The upper plot shows the up-spin and down-spin Mn 3d densities of states above and below the x-axis respectively, 
and the lower plot shows the up- and down-spin oxygen densities of states. In this case we have plotted only one of 
the t2g and one of the Cg bands for clarity. The only point that we want to make here is that the densities of states 
are very similar to those of the antiferromagnetic structure. Therefore as long as spin polarization is allowed, the 
electronic structure is not very sensitive to the relative polarization of neighboring Mn ions. 

b. Hexagonal The paraelectric hexagonal structure contains two Mn ions per unit cell, allowing studies of FM 
and AFM spin orientations without employing supercells. We find that the AFM phase is slightly lower in energy 
than the FM phase, as is the case in experiment. We focus on the AFM phase in what follows although, as we saw 
already in the cubic case, the spin ordering is not essential for our arguments. The properties that we consider here 
are determined by the local (i.e. internal to the oxygen cage) spin polarization rather than by long-range magnetic 
ordering. 

In contrast to the non-spin-polarizcd phase, in which the dxz, dyz and dxy, dx2-y2 doublets were roughly half 
occupied, here the Stoner exchange produces a ~ 2.5 eV energy splitting that causes these four orbitals to become 
almost completely spin polarized (Figure M). Their total up and down charges are indeed 3.93 and 0.47 electrons, 
respectively, while the occupation numbers of d 2 and d 2 orbitals are 0.53 and 0.30 electrons. The resulting magnetic 
moment on Mn (3.7 /is) is less than the 4 /is Hund's rule value because of the hybridization of d^ with Ot Pz 
states, and of dxy and dx2^y2 with Op Px and py states. However it is larger than the value calculated for the cubic 
case because of the larger exchange splitting occurring here. A non-negligible magnetic moment, M=0.15 fiB, is also 
present on Op, driven by exchange with the dxy and dx2-y2 orbitals, whereas the magnetization of Ot is almost 
vanishing, due to the AFM symmetry (for u—0 it would be exactly zero). 

The large spin splitting does not manage to open a gap. Indeed there is a tiny DOS at Ep, to which dxy and 
dx2-y2 orbitals from Mn, and Px and Py orbitals from Op contribute. This means that the bands crossing Ep come 
completely from orbitals oriented in the x ~ y plane (the tilting of the oxygen cages in the ferroelectric structure is 



probably sufficient to open an energy gap) . The mainly planar distribution of the charge density is also evident from 
the band energies shown in Figure |8|, calculated within the kz=Q plane of the Brillouin zone (left panel), and along 
[0001] (right panel). Only one band crosses Ep at kz—0, and its orbital character is a mix of d], and dl,2_ 2 states 
from Mn, and p]. and pj states from Op. The band energies are extremely flat along [0001], as is typical for a strongly 
layered compound. The system can be described as a poor metal in the x — y plane, and as an insulator perpendicular 
to it. If we discard the first unoccupied bands coming from in-plane orbitals, we obtain a one-dimensional 'pseudogap' 
of '^ 2 eV between occupied Ot P and empty Mn d^^ states, similar to the energy gap present between occupied 
oxygen states and empty transition metal d states in classic ferroelectric perovskite oxides such as BaTiOa. 

3. Effects of ferroelectric distortion on the electronic properties of cubic and hexagonal structures 

Finally we investigate the sensitivity of the charge distribution to atomic displacements by moving the Mn and 
in-plane oxygen atoms in opposite directions along the c and z axes in hexagonal and cubic phases respectively. 

The results for the hexagonal structure are shown in Figure S. Note that since the electric polarization occurs along 
the c axis, and the bipyramidal tilting docs not alter significantly the features of this compound, the paraelectric 
structure without the addition tilting is a sufficient starting point for the purpose of studying the nature of the 
ferroelectric displacement. For clarity, only the Mn orbitals most involved in the change of hybridization (i.e. d^z 
and dz2) are shown. The change of DOS upon distortion is very evident. The Mn d DOS is shifted down by ^ 2 
eV and is strongly hybridized with the Ot P states that are localized in a ^1 eV wide energy region. The atomic 
displacements increase the overlap of 0^22 and Ot Pz orbitals, inducing a pda- hybridization, and the overlap of d^z 
and dyz with Ot Px and Py orbitals, increasing the prf^ hybridization. Note that very similar changes in hybridization 
with ferroelectric displacement are seen in the prototypical ferroelectric BaTiOa [^. However in BaTiOa all the Ti 
3d orbitals are formally empty, and so displacement in any direction causes an energy lowering rehybridization as the 
oxygen ions transfer charge into the empty d orbitals. By contrast, the ferroelectricity in YMnOa is strongly uniaxial 
because only the 0^22 orbital is empty. Since the majority spin dxz and dyz orbitals are almost completely filled, their 
changes of hybridizations with oxygen on displacement, although significant, do not contribute to a change in electric 
polarization. 



The markedly different results which we obtain for the cubic structure are shown in Figure 10, Again only those 
orbitals most affected by the distortion (Mn ^^2 and d^z and Ot) are shown. The changes in DOS on distortion are 
much less pronounced and quite different in character from those which we found for the hexagonal phase. The most 
obvious change is a destabilization of the dxz band which results from the larger crystal field repulsion as the Mn is 
moved closer to the transverse oxygens along the z axis. Since the d^z orbitals are part of the filled i2g manifold they 
are unable to lower their energy by rehybridization with the oxygen 2p states, therefore the only contribution to their 
energetics is the unfavorable crystal field repulsion. Note that a typical Jahn- Teller distortion, in which the oxygens 
along the z axis move away from the cation, would have the opposite effect. In addition there is a slight broadening 
and shift down in energy of the d^-i band, which, since it was originally partly unoccupied, is now better able to accept 



charge density from the oxygen hgand that it has moved closer to. There is an almost rigid shift of the oxygen 2p 
bands to lower energy relative to the Fermi energy. This in fact a consequence of the shift of the Fermi level to higher 
energy due to the destabilization of the Mn 3d bands. The fact that the form of the oxygen bands does not change 
significantly is an indication that rehybridization is minimal. 

B. Comparison of the electronic properties of cubic YMnOa and LaMnOs 

Finally we briefly compare our calculated band structures of YMnOs in its metastable cubic perovskite structure, 
with those of the prototypical cubic manganite, LaMnOs. The purpose of this section is to determine whether there 
are any important differences between the electronic properties of cubic YMnOs and LaMnOs which might explain 
why the cubic phase is stable in LaMnOa but metastable in YMnOs. (Note that the cubic phase does not support 
ferroelectricity for either compound, therefore since LaMnOs does not form in the hexagonal structure it is never 
ferroelectric.) In fact we will find that the band structures of cubic YMnOa and LaMnOs are very similar, and 
therefore conclude that LaMnOs only forms in the cubic phase, whereas YMnOa can be either cubic or hexagonal 
(and therefore ferroelectric) because the yttrium ion is smaller than the lanthanum ion. 

First we report in Table || the relative energies of the paramagnetic (PM) , ferromagnetic (FM) and A-type ( AFM) 
phases of cubic YMnOs and LaMnOs. The relative energies of FM and AFM phases of LaMnOs are taken from 
Ref. |T2] and those for YMnOs were calculated in this work. The most important feature for our purposes is that 
the relative orderings of the states is the same in YMnOs as in LaMnOs. It is also striking that in both cases the 
fcrromagnctically ordered state is the most stable phase. This has been pointed out previously for LaMnOs 111], and 
the relationship between the observed Jahn- Teller distortion and the resulting A-type antiferromagnetic ground state 
is well understood |^l|]. Also, in both cases the calculated magnetic moment is about 3.0 /ib per Mn for the cubic FM 
and AFM structures, lower than the spin-only value of 4 fiB for the Mn'^^ ion as discussed earlier. The similarity in 
magnetic moment indicates similar amounts of spin-polarization and hybridization in both cases. 

Figure ^ shows the band structures of cubic PM YMnOs and LaMnOs plotted along the high symmetry axes of 
the simple cubic Brillouin Zone. The Fermi level is set to eV in both cases. The broad O 2p bands between around 
-2 and -7 eV can be seen clearly in both materials, with the Mn 3d bands above them, separated by an energy gap. In 
LaMnOs, the Mn 3d bands are largely separated in energy at each fc-point from the next highest energy bands, which 
are the La 5d bands. The highlighted lines in the LaMnOs band structure plots accentuate the upper and lower Mn 
3d bands, which have a similar form to each other and have been seen in other cubic perovskite manganites (see Refs. 
||ll| and |12|] .) This indicates a 'universality' in the manganite structure throughout the cubic perovskite manganite 
series, which is independent of the identity of the large cation. The same structure for the upper and lower Mn 3d 
bands can also be identified in the YMnOs band structure, however the Mn 3d bands are slightly broader than those 
in LaMnOs because the smaller lattice constant results in greater Mn 3d - O 2p overlap. The smaller lattice constant, 
and correspondingly larger crystal field repulsion, is also responsible for the larger t2g - eg splitting which can be 
clearly seen at the F point. In addition the Y 4d bands are lower in energy than the La 5d bands. These factors result 



in the lower part of the Y Ad bands intersecting the upper part of the Mn 3d bands in some regions of the Brillouin 
zone, making the universal structure less well-defined. 

To quantify any differences between cubic PM YMnOa and LaMnOs, we performed tight- binding analyses of the T 
to X regions of the respective band structures. Tight- binding parameters were obtained by non- linear- least-squares 
fitting |40|] to the calculated ab initio energies at the high symmetry F and X points, and at 19 points along the A 
axis. The tight-binding parameters thus obtained are given in Table 0. In both cases a good (root mean square (RMS) 
deviation from the ab initio energies of 0.20) fit was obtained by including only oxygen 2s and 2p and Mn 3d orbitals 
in the basis set, since we find no significant covalent bonding between the A cation and either the oxygen 2p or the 
Mn 3d orbitals. This is consistent with an early proposal by Goodenough M] that the magnetic properties of the rare 
earth manganites are determined by the Mn 3d - O 2p hybridization only. The parameters for YMnOa are similar 
to those for LaMnOa, except for an increase of around 10 % in the Mn 3d - O nearest neighbor overlaps. This is a 
consequence of the smaller lattice constant in YMnOa, and results in the slightly broader Mn 3d bands noted above. 

To test the effect of size on the electronic properties, we also calculated the properties of cubic paramagnetic YMnOa 
at the experimental LaMnOa lattice constant of 3.95 A . The resulting band structure of this hypothetical structure 
(not shown) is qualitatively similar to that of YMnOa at the true experimental lattice constant, except for a slight 
narrowing of the bandwidths as a result of the reduced overlap between the atomic wavef unctions, and a reduced t2g 
- Cg splitting as a result of the lower crystal field repulsion. 

We find that introduction of spin polarization causes the same kinds of changes in both compounds, therefore the 
conclusions drawn from analysis of the paramagnetic band structures continue to be valid. Figure O shows the up 
and down-spin band structures for YMnOs and LaMnOa along the high symmetry axes of the simple cubic Brillouin 
Zone. (We chose the ferromagnetic (FM) phases for simplicity in interpreting the band structures; as discussed above, 
we are interested in the effects of spin polarization rather than the nature of the ordering between magnetic moments.) 

First we examine the up- and down-spin LaMnOa band structures, and compare with the PM LaMnOa band 
structures to determine the changes which spin polarization causes in a 'conventional' manganite. The states which 
correspond to non-magnetic atoms arc unchanged from the paramagnetic state, and are identical for up- and down- 
spin electrons. For example, the dispersion of the lowest O 2p band is very similar for up- and down-spin, and for 
the PM phase. Also the La 5d bands, which were above the Mn 3d bands in the PM state are unchanged in form 
and energy. The characteristic perovskite manganite Mn 3d pattern which we remarked on earlier persists in the 
FM phase, appearing around 2 eV higher for the down-spin electrons than for the up-spin because of the exchange 
splitting. The up-spin Mn 3d and O 2p bands are strongly hybridized and there is no gap between them. However, 
the down-spin Mn 3d are split off from the O 2p bands by a larger gap than in the PM case. As a result, the Mn 3d 
bands occupy the same energy region as the unoccupied La 5d bands. The Fermi level cuts through the very bottom 
of the down-spin Mn 3d bands, and the conduction band is occupied almost entirely by up-spin electrons, indicating 
half-metallic behavior. 

Finally we compare the up- and down-spin YMnOa band structures. The first notable feature is that the non- 
magnetic Y Ad bands and O 2p bands are very similar for both up- and down- spin. This is because the exchange 



interaction in this hypothetical ferromagnetic YMnOs structure does not remove the gap between the majority Mn d 
bands and the O 2p band. (By contrast, in LaMnOs the higher O 2p bands overlap with the majority Mn 3d bands 
and therefore show a polarization dependence). This is again the result of the smaller lattice constant in YMnOs, 
which destabilizes the unfavorable FM interactions compared with those in LaMnOa, where the Mn ions are further 
apart. Both up- and down-spin Mn Sd bands show the characteristic perovskite manganite features (most noticeable 
in this case along the M-R-X directions) which we noted earlier in our discussions of the PM band structures. 

The strong similarities between the calculated properties of cubic YMnOa and LaMnOs that we have presented 
in this section confirm that the observed differences between the properties of hexagonal and cubic manganites are 
merely the result of the different sizes of the large cations in the two materials stabilizing different structures, and 
not due to any additional covalent bonding effects. 

V. CONCLUSIONS 

In summary, the formal similarity of d^ Mn'^+ in hexagonal and cubic YMnOa does not prevent the two phases from 
showing quite different properties. In particular, the difference between hexagonal and cubic field symmetry induces a 
very different charge distribution. In cubic pcrovskites, the Mn^+ ion has partially occupied e^ states. This creates an 
instability relieved by Jahn- Teller distortion and the ordering of Mn orbitals on the (001) plane. Both Jahn- Teller and 
orbital ordering tend to favor structural configurations that keep Mn and O aligned on the same plane, and therefore 
do not result in fcrroelectricity In the hexagonal perovskites, the crystal field produces an ordering of d states that 
leaves the dz2 orbital mostly unoccupied, and thus able to hybridize with the pz orbital of the transverse oxygen. This 
charge environment favorable to electric polarization along the c axis is only realized in condition of spin-polarization, 
not for the non-magnetic hexagonal phase. Thus, in hexagonal YMnOs the spin-polarization is actually the factor 
which enables the ferroelectricity. Finally, we emphasize that the transition from the AFM to the paramagnetic phase 
above T^r^SO K (i.e. well below the critical temperature of the ferroelectric phase), does not invalidate our picture, 
since it is the local spin-polarization of Mn'^+ that matters, independently of the actual presence (or absence) of spin 
ordering. Note also that there are no fundamental differences between the electronic structures of mctastable cubic 
YMnOa and the prototypical stable cubic manganite, LaMnOa. Therefore the metastability of the cubic perovskite 
phase in YMnOa is simply the result of the smaller size of the A cation which is unable to stabilize the cubic phase. 
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Figure Captions 

FIG. 1. The ideal cubic perovskite structure. The small B cation (in black) is at the center of an octahedron of oxygen 
anions (in gray). The large A cations (white) occupy the unit cell corners. 

FIG. 2. The ideal hexagonal perovskite structure. The small B cation (in black) is five-fold coordinated by oxygen anions 
(in gray). The dashed lines indicate the oxygen pentahedra. The large A cations (white) occupy seven- fold coordinated sites. 

FIG. 3. Calculated density of states for cubic paramagnetic YMnOs. 

FIG. 4. Orbital resolved densities of states for Mn and transverse oxygen in the paraelectric non-magnetic phase of hexagonal 
YMnOa. Only one DOS for each of the nearly degenerate pairs of orbitals {dxy,d^2_y2) and (dxz, dyz) is shown. 

FIG. 5. Orbital resolved densities of states on single Mn and O atoms in cubic antiferromagnetic YMnOa. 

FIG. 6. Orbital resolved density of states in cubic ferromagnetic YMnOs. 

FIG. 7. Orbital resolved DOS of a single Mn ion (top panel), in-plane oxygen Op (middle), and transverse oxygen Ot 
(bottom) in the paraelectric AFM phase of hexagonal YMnOa. p" and p'' indicate up and down components of O p DOS, 
respectively. 

FIG. 8. Band energies of paraelectric AFM hexagonal YMnOs. Q=7r/a (0,2/3^/^,0), P=7r/a (l,l/3^''^0)=^. 

FIG. 9. Orbital-resolved DOS for single Mn and O in hexagonal AFM YMnOs. The upper panel is for the paraelectric 
structure, the lower for the structure with Mn and O displaced along the c axis, in the fashion of a ferroelectric distortion. 

FIG. 10. Orbital-resolved DOS for single Mn and O in cubic A-type AFM YMnOs. The upper panel is for the paraelectric 
structure, the lower for the structure with Mn and O displaced along the z axis, in the fashion of a ferroelectric distortion. 

FIG. 11. Calculated band structures for cubic paramagnetic LaMnOa and YMnOs along the high symmetry axes of the 
Brillouin Zone. The highlighted lines in the band structure plots accentuate the upper and lower Mn 3d bands, which have a 
similar form to each other and to those of other perovskite manganites. 

FIG. 12. Up- and down-spin band structures for cubic LaMnOs and YMnOs along the high symmetry axes of the Brillouin 
Zone. 





LaMnOs 


YMnOa 


Eo2s 


-17.836376 


-17.9649 


Eo2p 


-4.514022 


-4.7382 


EMn-id 


-1.191285 


-1.4097 


Vo2s-02s 


-0.243683 


-0.2461 


V(02p~02p)iT 


0.620568 


-0.7167 


V(02p-02p)-r: 


-0.063466 


-0.1069 


\V{02p-02p)A'2 


0.182635 


0.1688 


[V{02p~02p)7t]2 


0.082951 


0.0620 


Vo2s~Mn3d 


-1.735814 


-1.9711 


V(02p~Mn3d)<y 


-1.838490 


-2.0601 


V(02p~Mn3d)Tv 


0.878961 


1.0837 


y(Mn3d-Mn3d)5 


0.066346 


0.0488 



TABLE I. Tight-binding parameters (in eV) for cubic LaMnOa and YMnOa obtained by non-linear-least-squares fitting to 
the ab initio eigenvalues along F to X. E indicates an orbital energy, and V an inter-atomic transfer integral. All transfer 
integrals are between nearest neighbors, except those with the subscript '2' which are between next nearest neighbors. Only 
the parameters listed in the table were allowed to be non-zero in the fitting procedure. 





YMnOs 


LaMnOs 


FM 

A-AFM 

PM 




+140 meV 
+0.70 eV 



+110 meV 
+ 1.16 eV 



TABLE II. Relative energies of different magnetic phases in cubic paramagnetic YMnOs and LaMnOs 
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